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A  USERS  GUIDE  TO  KLYSMA  II:  A  MULTIFLUID  TRANSPORT  CODE 


I.  Introduction 

This  report  documents  a  new  one  dimensional  multlfluld  transport  code, 
KLYSMA  II,  designed  to  study  the  cross-field  transport  of  magnetized  plasma  in 
the  presence  of  both  anomalous  and  classical  transport  processes,  radiation, 
and  chemical  reactions*  Because  the  code  was  designed  for  research,  as  well 
as  production  purposes,  versatility  was  incorporated  at  its  conception.  It  is 
modular,  well  documented,  and  incorporates  a  variety  of  logical  switches  so 
that  many  physical  processes  may  be  switched  on  and  off  to  study  their  effect 
on  the  problem  under  consideration.  In  addition  the  code,  while  built  on  a 
VAX-780,  was  constructed  with  vectorization  in  mind.  We  note  that  KLYSMA  II 
differs  from  KLYSMA  in  many  ways:  better  numerical  algorithms  are  used  (e.g., 
advanced  FCT,  implicitly  solves  the  coupled  set  of  ion  and  electron  diffusion 
equations),  new  physics  is  Incorporated  (e.g.,  electrons  and  ions  are 
thermally  coupled,  new  instabilities  are  considered  such  as  the  lower-hybrid- 
drift  instability),  the  code  is  vectorized,  and  the  code  is  well-documented. 
In  addition,  very  little  physics  is  "hard  wired"  in  KLYSMA  II  so  that 
modifications  and/or  upgrades  can  be  easily  made  as  the  need  arises. 

In  section  II  we  present  the  physical  processes  Incorporated  into  KLYSMA 

II,  the  equations  that  describe  these  processes,  the  transport  coefficients 
used,  and  the  numerical  algorithms  utilized.  In  section  III  we  present  the 
structure  of  KLYSMA  II  and  describe  the  purpose  of  each  module.  In  section  IV 
all  the  logical  control  parameters  are  defined  as  well  as  their  purpose. 
Further  a  definition  for  each  global  array  or  parameter  utilized  by  KLYSMA  II 
will  be  given.  Finally  in  section  V  we  present  a  model  run  to  illustrate  the 
use  of  the  code. 

II.  Basic  Physical  Processes  and  Numerical  Algorithms 

KLYSMA  II  Is  designed  to  study  how  the  cross-magnetic  field  transport  of 
a  species  of  ions  j  with  density  nj,  temperature  Tj  and  velocity  Vj  Interacts 
one  dimensionally  with  an  ambient  plasma  of  ion  density  n^,  ion  temperature 
T^.and  velocity  V^,  electron  density  nfi,  electron  temperature  Te,  and  electron 
velocity  Vg  in  the  presence  of  a  magnetic  field  _B.  Because  the  heating  and 
momentum  coupling  rates  may  be  controlled  by  either  collective  ("anomalous") 
plasma  processes  and/or  classical  binary  collisions,  prescriptions  for  the 

computation  of  collective  cross-field  transport  processes,  as  well  as 
Manuscript  approved  March  2,  1984. 
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classical  transport  processes,  are  incorporated  into  KLYSMA  II.  The  anomalous 
transport  coefficients  and  the  prescriptions  for  their  use  that  are 
Incorporated  into  KLYSMA  II,  with  the  exception  of  the  lower  hybrid  drift 
instability  (LHDI)  and  the  drift  cyclotron  instability  (DCI),  can  be  found  in 
Lampe  et  al  (1975).  The  transport  coefficients  resulting  from  the  LHDI  and 
DCI,  and  the  prescriptions  for  their  use  were  obtained  from  Huba  (1983).  All 
classical  transport  coefficients  utilized  in  KLYSMA  II  can  be  found  in 
Braglnskil  (1965). 

The  following  equations  are  solved: 

3n . 

3^  +  1  .(njVj)  -  0  ,  (1) 

ion  continuity  equation, 


it  <VA>  +  2-<*jnAV  * 

-  +  ziV<E  +  +  • 


(2) 


ion  momentum  equation, 
3P 
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ion  pressure  equation. 
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(4) 


force  balance  on  electrons. 
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electron  pressure  equation, 
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quasi-neutrality 


Faraday's  Equation, 
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Ampere's  Equation, 


(9) 


where  Pj  ■*  Ojk^Tj,  e  is  the  magnitude  of  the  electron  charge,  Zj  the  charge 
state  of  the  ion  species  j,  mj  the  mass  of  ion  species  j,  and  me  the  electron 
mass. 

Rj  represents  the  rate  of  momentum  exchange  between  ion  species  k  and  j, 
and  between  ion  species  j  and  electrons,  and  has  the  general  form 


nkmknjmj 


R  m  V  v  K  K  J  J  (V  -  V  ) 

-J  Li.e  jk  njmj+nkmk  -4t  “j 


(10) 


where  Vj^  represents  the  total  effective  collision  frequency  (both  anomalous 
and  classical)  between  species  k  and  ion  species  j  or  the  electrons  .  Qj 
represents  the  rate  at  which  ion  species  j  is  heated  or  cooled  cm-^s-^  due  to 
energy  exchange  between  ion  species  j  and  the  other  plasma  constituents  and 
has  the  form 


_  n,m  ,n.  m. 

Q,  -  l  v,.  (V,  -  V.)2  j  j  k  k -  +  — - 

j  t  jk  -*  -j  njBj  +  n^  Bj 


3m  n  k. 

*  -V1  <Te  '  V- 
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where  the  index  k  in  the  sum  runs  only  over  the  ion  species, 

3/m  (k  T  )3/2 
e  b  e 

Tej  "  4/2tt  Xe^Z^nj  ’ 


(ID 


and  X  is  the  Coulomb  logarithm  (Braginskii,  1965,  p.215).  represents  the 
perpendicular  to  J5  heat  flux  carried  by  the  ion  species  j  and  has  the  form 


Y  ^  O 


T  > 


j  b  j  1  (ax2  +fl) 


x*  +  yx^  +  5 


3/m  (k.  T  )3/2 

-  MINIMUM  t-ttnhk—' ■  T  jan^  * 


Tjan  being  the  anomalous  value  of  Tj  due  to  collective  processes,  x  ■ 
ncJTj,  -  Z^eB/mjC  ,  and  a,  8*  y,  6  are  given  by  Braginskll 

according  to  the  charge  state  of  the  Ion  species  j. 
qe  Is  the  electron  heat  flux  given  by 


(1965) 


qe  -  -KJ  VjdCh^), 


where 


a«VeTeff  <V  +  V 

mg  (y4  +  djX2  +  6q)  * 


x  m  MINIMUM  (t  t  ),  x  being  the  anomalous  value  of  t  ,  y  ■ 
etr  e,  an  an 

nceTeff,  nce  "  eB/mec*  and  Yjt  Yc»  «c  are  given  by  Braginskll. 

The  quantity  Rg  and  Qg  represent,  respectively,  the  rate  at  which 
momentum  and  energy  Is  gained  or  lost  by  electrons  and  have  the  form 


-  I  Weff'Ik  -  V 


3m  n  k, 

Q  -  l  {m  n  v  ,C(V.  -  V  )2  +  — 6  e  (T,  -  T  )}  -  R  ., 
■^e  £  1  e  e  effv-k  -e'  m^  Tgk  v  k  e/J  rad’ 


where  Rra(j  represents  radiation  losses. 

These  equations  are  solved  In  a  time  split  fashion  because  of  the  wide 
range  of  time  scales  associated  with  the  physical  processes  contained  In 
equations  (I)  -  (9).  Momentum  exchange  between  Ion  species,  and  electrons 
represented  by  jtj  In  (2),  is  treated  implicitly,  as  are  the  thermal  transport 
and  the  source  and  sink  terms  generically  denoted  by 


in  equations  (3)  and  (S).  All  other  source  terms  in  equations  (2)-(5)  are 
treated  as  explicit  source  terms  for  the  appropriate  calls  to  the  explicit  FCT 
module  DSSFCT  (Boris,  1976). 

The  module  MULTIF  performs  all  convective  transport  by  splitting  the 
Integration  into  half  and  full  time  steps  so  as  to  achieve  space  and  time 
centering.  Module  MULTIF  calls  MONTRN  which  implicitly  exchanges  momentum 
among  the  various  plasma  constituents.  If  the  user  selects  LTCENT  -  .FALSE, 
the  hydro  tlmestep  DT  is  split  into  many  subcycles  each  with  a  time  step  DTC 
such  that  DT  *  NTIME  *  DTC.  DTC  is  the  time  step  associated  with  the  shortest 
physical  process  contained  in  Q  of  (16).  Under  these  circumstances  the 
classical  and  anomalous  transport  coefficients  are  updated  NTIME  times  with 
DIFFX3  being  called  after  each  update.  DIFFX3  is  an  implicit  module  designed 
to  be  fully  compatible  with  the  FCT  module.  Physically  subcycling  the 
heating,  cooling,  and  equilibration  terms  in  the  pressure  equations  is  fully 
justified  because  Q  changes  very  rapidly  on  a  hydro  time  scale  making  it 
necessary  to  update  the  flow  quantities  less  frequently.  However,  in  the 
event  the  user  prefers  a  fully  time  centered  run,  which  is  less  economical  but 
nevertheless  more  accurate,  taking  LTCENT  -  .TRUE,  will  achieve  this  so  that 
DT  -  DTC  and  NTIME  -  1.  DIFFU3,  which  is  called  by  DIFFX3,  solves  (implicitly 
the  ion  and  electron)  diffusion  equations,  together  with  equilibration  and 
source  terms,  simultaneously  utilizing  the  time  step  determined  by  DTSET  and 
the  condition  of  LTCENT.  It  incorporates  a  special  finite  differencing 
technique  that  requires  the  energy  flux  from  the  center  of  the  cell  I  to  the 
interface  of  the  cell  to  equal  the  energy  flux  from  the  interface  to  the 
center  of  the  cell  1+1.  This  allows  the  temperature  drop  between  cells  to 
occur  predominately  in  the  cell  with  the  smaller  coefficient  of  thermal 
conduction,  as  would  occur  physically. 

Subroutine  REZONE  allows  the  user  to  "window"  the  flow  via  a  sliding 
rezone  technique  that  uses  the  maximum  and  minimum  velocity  of  ion  specie  2 
(the  debris  or  "injected  Ion”)  to  slide  the  window  along  with  the  flow  if  FLAG 
*  0.0,  or  a  Lagrangian  capability  if  FLAG  ■  1.0.  The  subroutine  DTSET 
computes  self-consistently  the  time-step  necessary  for  numerical  stabiity  as 


well  as  enforcing  the  condition  that  no  physical  quantity  vary  by  more 
than  ~  5%  In  a  timestep. 

Ideally,  debris  (Ion  fluid  2)  -  air  (ion  fluid  1)  coupling  should  be 
studied  with  a  macroparticle  or  hybrid  scheme  which  resolves  the  complex 
interaction  of  the  two  streams  in  velocity  space.  This  is  not  always  feasible 
for  reasons  which  Include  economics  and  machine  limitations.  A  compromise  can 
be  achieved  through  the  use  of  a  multi-fluid  scheme.  It  permits  the  air  and 
debris  plasmas  to  maintain  separate  identities,  and  to  interpenetrate.  With  a 
single  fluid  scheme,  stagnation  flow  occurs  at  the  debris-air  interface  and 
all  local  relative  streaming  energy  is  transformed  into  thermal  energy.  There 
is  no  way  to  study  anomalous  ion-ion  coupling  with  the  latter  scheme,  since 
the  single  fluid  assumption  Implies  a  single  fluid  velocity  at  every  point. 

With  as  few  as  two  fluids,  a  multifluid  scheme  can  reproduce  rather 
faithfully  the  configuration  in  phase  space  predicted  by  a  macroparticle 
scheme.  In  one  respect,  however,  it  will  eventually  predict  incorrect 
behaviour  unless  corrective  measures  are  taken.  Air  plasma  which  has 
"coupled"  via  collisions,  classical  or  collective,  or  via  Larmor  processes,  to 
the  expanding  debris  plasma  subsequently  moves  along  with  the  debris,  i.e.,  it 
is  picked  up.  This  picked  up  air  and  debris  mixture  then  interacts  with  the 
incoming  upstream  air  plasma,  via  the  same  processes  which  earlier  coupled  the 
debris  to  the  picked  air,  rather  than  stagnating  and  producing  a  shock.  One 
way  to  deal  with  this  pickup  is  to  reassign  some  portions  of  the  picked-up  air 
fluid  to  the  debris  fluid.  This  is  implemented  in  KLYSMA  via  SHUFF.  SHUFF 
uses  the  assumption  that  both  air  and  debris  have  Maxwellian  radial  velocity 
distributions.  It  then  imposes  a  velocity  cutoff,  such  that  all  plasma  above 
the  cutoff  velocity  is  considered  to  be  debris,  and  all  plasma  below  It  Is 
taken  to  be  air.  In  addition,  conservation  of  mass,  momentum,  energy,  and 
charge  are  imposed  in  the  transfer  process. 
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HI.  Code  Structure 


As  noted  in  the  introduction  the  code  is  modular,  the  code  is  driven  from 
the  main  subroutine  MAINMF,  which  is  structured  as  follows: 


-  initiates  all  arrays 
and  boundary  conditions 


-  restarts  code  as  required 


—  Computes  necessary 

diagnostics  and  outputs; 
plots  and  tabulates 
data 


-  Chemistry  is  advanced 


-  Ion  leak  losses  are 
computed 


—  computes  necessary 
parameters  for  use  in 
ANOMAL  and  CLASS 


-  computes  anomalous 
transport  coefficients 


STOP 

END 


-  computes  classical 
transport  coefficients 


-  Equilibrates  ions  and 
electrons  and  diffuses 
ions  and  electron  heat 
fluxes  together  with 
thermal  sources  and  sinks 
implicitly  using  DTC 
obtained  from  DTSET 


-  Performs  all  convective 
transport  explicitly  using 
DT  obtained  from  DTSET 


-  sets  to  zero  all 
necessary  arrays 


-  Computes  time  step  as 
necessary  for  numerical 
stability  both  DT  and  DTC 
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SUBROUTINE  INITIAL 


This  subroutine  reads  initial  data,  boundary  conditions,  and  boundary 
condition  factors  off  disc  to  initiate  the  run.  In  addition,  it  sets  up  the 
initial  conditions.  In  general  this  fortran  nodule  will  be  the  one  most 
likely  modified  by  the  user- 

SUBROUTINE  RSTART 

This  subroutine  restarts  the  run  at  the  necessary  time  step  by  reading 
from  disc  the  data  written  there  by  DUMPER.  Subroutine  DUMPER  must  be  called 
at  least  once  before  calling  RSTART.  For  example,  the  user  may  want  to  run 
the  code  first  from  1  to  the  2001  timestep,  and  then  dump  the  data  to  disc  for 
restarting.  Thus  initially  LDUMP  =  .TRUE,  and  LRESTR  =*  .FALSE.  When 
restarting  the  code  at  timestep  2001  LRESTR  must  be  .TRUE..  LDUMP  may  be 
•TRUE,  or  .FALSE,  depending  on  whether  the  user  wants  to  restart  KLYSMA  II 
again. 

SUBROUTINE  ZERO 

This  subroutine  simply  sets  to  zero  all  arrays  that  are  used  in  summation 
processes. 

SUBROUTINE  COUPLE 

This  subroutine  computes  a  variety  of  physical  parameters  needed  to 
compute  both  the  classical  and  anomalous  transport  coefficients.  Both  ANOMAL 
and  CLASS  are  entries  to  COUPLE.  COUPLE  must  be  called  prior  to  either  ANOMAL 
or  CLASS. 

ENTRY  ANOMAL 

This  entry  into  couple  computes  the  cross-field  anomalous  transport 

coefficients  based  on  Lampe  et  al  (1975). 

ENTRY  CLASS 

This  entry  into  couple  computes  the  cross-field  classical  transport 

coefficients  obtained  from  Braginsky  (1965). 


SUBROUTINE  PRNT 


This  subroutine  prints  all  transported  quantities  as  well  as  the 
individual  source  terms  for  the  momentum  and  temperature  equations.  I PRNT 
determines  the  number  of  time  cycles  between  calls  to  PRNT. 

SUBROUTINE  PAGE PL  (CALCM) 

This  subroutine  generates  page  plots  (versatex  plots  using  NCAR  plot 
package)  of  the  same  quantities  printed  by  PRNT.  IPAGE  determines  the  number 
of  time  cycles  between  calls  to  PAGEPL  (CALCM). 

SUBROUTINE  DIAG 

This  subroutine  generates  a  time  (or  ordinate)  versus  distance  (abscissa) 
curve  of  the  position  of  the  coupling  shell.  The  point  at  which  the  magnetic 
field  has  a  maximum  is  used  as  the  spatial  position  of  the  coupling  shell  at 
time  t. 

SUBROUTINE  DTSET 

This  subroutine  is  designed  to  determine  the  maximum  timestep  permitted 
that  will  still  be  numerically  stable.  It  computes  both  DTC  and  DT. 

SUBROUTINE  DIFFX3 

This  subroutine  transports  the  thermal  heat  flux  of  the  electrons  and 
ions  perpendicular  to  the  attendant  magnetic  field  as  well  as  the  thermal 
sources  and  sinks.  It  utilizes  a  special  subroutine  DIFFU3  that  implicitly 
solves  the  coupled  set  of  ion  and  electron  diffusion  equations  and  is  designed 
to  be  fully  compatible  with  the  FCT  Module  DSSFCT.  A  call  to  DIFGE3  must  be 
made  prior  DIFFU3  so  that  the  updated  metric  coefficients  needed  in  DIFFU3  are 
generated . 

SUBROUTINE  MULTIF 

This  subroutine  is  the  heart  of  KLYSMA.  It  utilizes  the  FCT  Module 
DSSFCT  to  solve  the  various  fluid  equations.  Note  that  the  transport 
coefficients  used  in  these  source  terms  are  not  time  centered  when  LTCENT  ■ 
•FALSE.  However,  when  LTCENT  ■  .TRUE,  all  of  the  transport  coefficients, 
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anomalous  and  classical,  are  recomputed  both  at  the  half  and  the  full  time 
step.  In  addition,  DIFFX3  Is  called  at  both  the  half  and  full  time  step  if 
LTCENT  »  .TRUE. 

MULTIF  also  calls  a  Subroutine  REZONE  that  keeps  the  flow  "windowed" 
spatially  between  the  maximum  and  minimum  values  of  the  debris  velocity. 

SUBROUTINE  SHUFF 

This  subroutine  is  designed  to  allow  air  (J=l  fluid)  to  be  picked  up  by 
the  faster  moving  debris  (J*2  fluid)  and  for  any  deaccelerating  debris  to  be 
converted  to  slower  air.  To  accomplish  this  an  algorithm  is  used  that  assumes 
both  ion  species  are  locally  Maxwellian  within  a  grid  cell  and  drifting  with  a 
speed  VIX(I,J).  In  addition,  it  is  assumed  that  whenever  the  air  fluid 
velocity  exceeds  VD  *  (VIX(I,1)  +  VIX(I,2))/2  that  fraction  of  the  air 
distribution  that  exceeds  VD  becomes  debris.  In  other  words  the  fraction 

vp 

J-oo  ^alr  dvx  stays  a*r  an<*  tlie  fracti°n  fv  ^air^vx  becomes  debris.  Likewise 
for  the  debris.  With  these  assumptions,  together  with  the  requirements  that 
mass,  momentum,  energy  and  charge  be  conserved  within  a  cell  volume  adjusted 
densities,  velocities,  temperatures,  and  charges  are  computed. 

SUBROUTINE  HANE$ 

This  subroutine  contains  three  entries,  HANE$,  SOLAR$,  and  LASER$ ,  each 
of  which  computes  the  initial  conditions  for  a  HANE,  SOLAR,  or  LASER  Run. 

SUBROUTINE  ILEAK 

This  subroutine  computes  the  total  number  of  debris  ions  loss  from  the 
coupling  shells,  as  well  as,  the  total  debris  energy  loss  based  on  the  loss 
cone  instability. 


IV.  DEFINITIONS 


The  Module  GLOBMF.FOR  contains  all  the  parameters  global  to  KLYSMA.  The 
following  common  blocks  are  included  in  GLOBMF.  FOR. 

COMMON/ DIM/ 

NSPEC:  Number  of  ion  species,  presently  dimensioned  as 

MSPEC  *  2  in  the  parameter  statement 

NCHEM:  Number  of  ionization  levels  followed  for  each  ion 

species 

NINSTE:  Number  of  electron-ion  two  stream  instabilities, 

presently  dimensioned  as  MINSTE  =  4  in  the  parameter 
statement 

NINSTI:  Numer  of  ion-ion  two  stream  instabilities,  presently 

dimensioned  as  MINSTI  =>  2  in  the  parameter  statement 


COMMON/ PHYSIC/ 

NI  (NPT,  MSPEC): 

TI  (NPT,  MSPEC): 
VIX  (NPT,  MSPEC): 
VIY  (NPT,  MSPEC): 
VIZ  (NPT,  MSPEC): 
PI  (NPT,  MSPEC): 
RHO  (NPT,  MSPEC): 
RHOI  (NPT,  MSPEC): 
NE  (NPT): 

VEX  (NPT): 

VEY  (NPT): 

VEZ  (NPT: 

BY  (NPT): 

BZ  (NPT): 

EX  (NPT): 

EY  (NPT): 

EZ  (NPT): 

TE  (NPT): 

JY  (NPT): 

JZ  (NPT): 

PE  (NPT): 


ion  number  densities  (particles  cm  ) 
ion  temperatures  (degrees  Kelvin) 
ion  drift  velocities  in  x  direction  (cm/s) 
ion  drift  velocities  in  y  direction  (cm/s) 
ion  drift  velocities  in  z  direction  (cm/s) 
ion  pressures  (dynes/cm  ) 

_o 

ion  mass  densities  (grams  cm  ) 

Inverse  of  RHO 

electron  number  density  (particles/cm  ) 
electron  drift  velocity  in  x-direction  (cm/s) 
electron  drift  velocity  in  y-direction  (cm/s) 
electron  drift  velocity  in  z-direction  (cm/s) 
Y-component  of  magnetic  field  (gauss) 
Z-component  of  magnetic  field  (gauss) 
X-component  of  electric  field  (statvolts/cm) 
Y-component  of  electric  field  (statvolts/cm) 
Z-component  of  electric  field  (statvolts/cm) 
electron  temperature  (degrees  Kelvin) 

2 

y  component  of  current  density  (statamps/cm  ) 

2 

z  component  of  current  density  (statamps/cm  ) 

_2 

electron  pressure  (dynes  cm  ) 
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UNIT  (NPT) : 
ZERO  (NPT): 
PEDX  (NPT): 


array  of  which  all  elements  are  one 
array  of  which  all  elements  are  zero 
derivative  of  PE  with  respect  to  x 


COMMON/ CYCLE 
MINSTP: 
MAXSTP: 
IDIAG: 
IPRNT: 

I PAGE: 
IFLICK: 
ISTEP: 
IDUMP: 


index  of  STARTING  VALUE  of  time  loop 
index  of  final  value  of  time  loop 
frequency  of  calls  to  DIAG 
frequency  of  calls  to  PRNT 
frequency  of  calls  to  PAGE PL 

frequency  of  calls  to  FLICK  (dummy  at  present) 
value  of  index  corresponding  to  timestep 
frequency  of  intermediate  dumps 


COMMON/MESH/ 
XCOR  (NPT): 
XCORH  (NPC) : 
NX: 

NX1: 

NX2 : 

ALPHA: 


FLAG: 


DT: 
TIME : 
DTC: 


NTIME : 
COURT: 


positions  of  cell  centers 
positions  of  cell  interfaces 
number  of  cells  (or  grid  points) 

NX+1  number  of  cell  interfaces 
NX+2 

Geometry  chosen  cartesian  (ALPHA  =  1),  cylindrical 

(ALPHA  =2),  or  spherical  (ALPHA  =  3) 

used  if  LAGRANGIAN  code,  FLAG  =0.0  Eulerian, 

FLAG  =  1.0  Lagrangian,  must  be  used  in  conjuction 
with  LREZONE  =  .TRUE, 
hydro  time  step  (secs) 

TIME  (secs) 

TIME  STEP  associated  with  shortest  physical  process 

modelled  by  code 

NTIME  =  MAX  (DT/DTC,  1) 

COURANT  FACTOR 


COMMON/ PARA/ 

CS  (NPT,  MSPEC) :  ion  sound  speed  (cm/s) 

WPI  (NPT,  MSPEC):  ion  plasma  frequency  (Hertz) 


13 


n 


WCI  (NPT,  MSPEC) : 
VTI  (NPT,  MSPEC): 
VA  (NPT,  MSPEC): 
VDEIX  (NPT,  MSPEC): 


VDEIP  (NPT,  MSPEC): 

RCI  (NPT,  MSPEC): 
GYRRAD  (NPT, MSPEC): 
UIJ  (NPT): 

VTE  (NPT): 

WPE  (NPT): 

WCE  (NPT): 


ion  cyclotron  frequency  (Hertz) 

ion  thermal  speed  (cm/s) 

individual  ion  Alfven  speed  (cm/s) 

absolute  value  of  ion-electron  drift  speeds  in 

x-direction  (cm/s) 

absolute  value  of  ion-electron  drift  speeds 

perpendicular  to  B_  (cm/s) 

ion  thermal  Larmor  radii 

ion  non-thermal  Larmor  radii 

absolute  value  of  ion-ion  drift  speed  in 

in  x-direction  (cm/s) 

electron  thermal  speed  (cm/s) 

electron  plasma  frequency  (Hertz) 

electron  cyclotron  frequency  (Hertz) 


COMMON/ BOUND/  see 
NIR  (MSPEC): 

NIL  (MSPEC): 
NIRO  (MSPEC): 
NILO  (MSPEC): 
VIXR  (MSPEC): 
VIXL  (MSPEC): 
VIXRO  (MSPEC): 
VIXLO  (MSPEC): 
VIYR  (MSPEC): 
UIYL  (MSPEC): 
VIYRO  (MSPEC): 
VIYLO  (MSPEC): 
TIR  (MSPEC): 

TIL  (MSPEC): 
TIRO  (MSPEC): 
TILO  (MSPEC): 


TERO: 

TELO: 


Boris  NRL  Memo  Report  (3237) 

RIGHT  ion  density  boundary  condition  factor 

LEFT  ion  density  boundary  condition  factor 

RIGHT  ion  density  boundary  value 

LEFT  ion  density  boundary  value 

RIGHT  ion  x-velocity  boundary  condition  factor 

LEFT  ion  x-velocity  boundary  condition  factor 

RIGHT  ion  x-velocity  boundary  value 

LEFT  ion  x-velocity  boundary  value 

RIGHT  ion  y-velocity  boundary  condition  factor 

LEFT  ion  y-velocity  boundary  condition  factor 

RIGHT  ion  y-velocity  boundary  value 

LEFT  ion  y-velocity  boundary  value 

RIGHT  ion  temperature  boundary  condition  factor 

LEFT  ion  temperature  boundary  condition  factor 

RIGHT  ion  temperature  boudary  value 

LEFT  ion  temperature  boundary  value 

RIGHT  electron  temperature  boundary  condition  factor 
LEFT  electron  temperature  boundary  condition  factor 
RIGHT  electron  temperature  boundary  condition 
LEFT  electron  temperature  boundary  condition 
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BZR: 


BZL: 

BZRO: 

BZLO: 


RIGHT  z-component  of  _B  boundary  condition  factor 
LEFT  z-component  of  B  boundary  condition  factor 
RIGHT  z  component  of  _B_  boundary  value 
LEFT  z-component  of  B  boundary  value 


COMMON/ CHEM/ 

RAD  (NPT):  Radiation  Losses  (ergs/cm^/s) 

NIZ  (NPT,  MSPEC,  NCHEM):  Particles  cm-^  for  each  ionization  state  of  a 
given  ion  specie 


COMMON/TRANS/ 

VPIX  (NPT,  MSPEC,  MINSTE) :  Phase  velocities  in  x-direction  for  a  given 
electron-ion  two  stream  instability  (cm/s) 

VPIEY  (NPT,  MSPEC,  MINSTE):  Phase  velocities  in  y-direction  for  a  given 
electron-ion  two  stream  instability  (cm/s) 

VEFF  (NPT,  MSPEC,  MINSTE):  effective  collision  frequencies  resulting 
from  a  given  electron-ion  two  stream  instability  (Hertz) 

VEFFII  (NPT,  MINSTI) :  effective  ion-ion  collision  frequency  resulting 
from  a  given  ion-ion  two  stream  instability  (Hertz) 

VPII  (NPT,  MINSTI):  phase  velocities  in  x-direction  for  a  given  ion-ion 
two  stream  instability  (cm/s) 

EQBIE  (NPT,  MSPEC):  classical  electron-ion  equilibration  frequencies 
(Hertz) 

VEI  (NPT,  MSPEC):  classical  electron-ion  collisional  frequencies  (Hertz) 
COLJE  (NPT,  MSPEC):  total  effective  electron-ion  collision  frequency, 
anomalous  +  classical  (Hertz) 

KAPI  (NPT,  MSPEC):  ion  thermal  conduction  coefficient 

VII  (NPT,  MSPEC):  classical  ion-ion  collision  frequency  (Hertz) 

EQBIK  (NPT,  MSPEC):  classical  ion-ion  equilibration  frequency  (Hertz) 

COLKJ  (NPT,  MSPEC):  classical  ion-ion  momentum  exchange  rate 
OFFLAR  (NPT,  MSPEC):  binary  switch  used  to  turn  off  Larmor  coupling 

terms  if  effective  ion-ion  collision  frequency  exceeds  ion  gyro  frequency 
COLJK  (NPT,  MSPEC,  MINSTI):  anomalous  ion-ion  momentum  exchange  rate 
RESIS  (NPT):  classical  resistivity  (sec) 

RAPE  (NPT):  electron  thermal  conduction  coefficient 
RSIZE:  system  size  parallel  to  B  field  (cm) 


COMMON/ TIMSTP/ 


DRAG(NPT,  16,  MSPEC) :  individual  terms  that  result  in  ion  acceleration  or 

drag  (dynes) 

IONH(NPT,  16,  MSPEC):  individual  terms  that  result  in  ion  cooling  or 

O 

heating  (ergs/cm  /s) 

ELEC1I(NPT,  16):  individual  terms  that  result  in  electron  heating  or 

cooling  (ergs/cm  /s) 


COMMON/CONST/ 

-  all  physical  constants  are  c.g.s. 

EC: 

Electron  charge,  ECI  =  1.0/EC 

BK: 

Boltzmann's  constant,  BKJ  =  1.0/BK 

C: 

Speed  of  light,  Cl  =  1.0/C 

PIE: 

3.145927 

PIE4: 

4.0  *  PIE 

PIE8: 

8.0  *  PIE 

EM: 

electron  mass 

EMEC: 

EM/EC 

GAME: 

electron  ratio  of  specific  heats 

GAMEl : 

GAME-1 

GAME 2: 

GAME-2 

GAMI 

ion  ratio  of  specific  heats 

GAMI1 : 

GAMI-1 

GAMI 2: 

GAMI-2 

EVK: 

Conversion  fractor  from  eV  to  degrees  Kelvin 

COMMON/ ION/ 

MASS  (MSPEC): 

ion  masses  (grams) 

MASSI  (MSPEC): 

inverse  of  MASS 

MR  (MSPEC): 

electron  to  ion  mass  ratio  for  each  ion  specie 

SQMR  (MSPEC): 

inverse  of  SQMR 

S MASS I  (MSPEC): 

square  root  of  MASSI 

Z(NPT,  MSPEC): 

ion  charge 

COMMON/LOGIC/ 


LPAGE : 

control 

switch  to 

call 

PAGE PL 

LPRNT: 

control 

switch  to 

call 

PRNT 

LCPLE : 

control 

switch  to 

call 

COUPLE 

LANOM: 

control 

switch  to 

call 

ANOMAL 

LCLASS : 

control 

switch  to 

call 

CLASS 

LRESTR: 

control 

switch  to 

call 

RSTART 

LDUMP: 

control 

swtich  to 

call 

DUMPER 

LNEU: 

control 

switch  to 

call 

NEUTRALS  (Dummy  for  now) 

LRAD: 

control 

switch  to 

call 

RADIAT  (Dummy  for  now) 

LDOAG: 

control 

switch  to 

call 

DIAG 

LIAI: 

switches 

on  ion-acoustic  instability  if  LANOM  = 

TRUE 

LBCI: 

switches 

on  beam  cyclotrm  instability  if  LANOM  : 

-  TRUE 

LMTI : 

switches 

on  modified  two  stream  instability  if 

LANOM  -  TRUE 

LMII: 

switches 

on  magnetized 

ion-ion  instability  if  LANOM  *  TRUE 

LUII : 

switches 

on  unmagnetized  ion-ion  instability  if 

LANOM  -  TRUE 

LHDI: 

switches 

on  unmagnetized  lower  hybrid  drift  instability 

if  LANOM  -  TRUE 

LVISCO: 

switches 

on  artificial 

viscosity 

LDIFF : 

control 

switch  for  DIFFX3 

LARMOR: 

switches 

;  on  or  off  ion  Larmor  coupling 

LSHUFF : 

control 

switch  to 

call 

SHUFF 

LILEAK: 

control 

switch  to 

call 

ILEAK 

LTCENT : 

if  true 

transport 

is  fully  time  and  space  centered 

LSTION: 

control 

switch  to 

call 

Supra  thermal  ion  module 

(dummy  for  now) 

LBOHM: 

control 

switch  to 

call 

BOHM  (dummy  for  now) 

LCHEM 

control 

switch  to 

call 

chemistry  module 

LJAC71 : 

control 

switch  to 

call 

JAC71  -  not  used 

LHANE 

control 

switch  to 

call 

HANE 

LSOLAR 

control 

switch  to 

call 

SOLAR 

LLASER 

control 

switch  to 

call 

LASER 

LECHO 

control 

switch  to 

call 

to  give  positions  for  module  transfer 

LSTOP 

control 

switch  to 

call 

to  stop  code 

LFLICK 

control 

switch  to 

call 

FLICK  (Dummy  for  now) 

LREZON 

control 

switch  to 

call 

REZONE 
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LCALCM 


control  switch  to  call  CALCM 
LPARTX  control  switch  to  call  PARTX  (Dummy  for  now) 

COMMON/ REZONE/ 

IMAX: 

IMIN: 

VNEW  (NPT) : 

XOLD  (NPT): 

COMMON/HANE/ 

HRA: 

LAT: 

DAY: 

ALT: 

YIELD: 

WMASS : 

RXRADN : 

BACKIZ: 

NTOTAL : 

RMASS: 


total  number  of  debris  ions  loss  down  field  lines 

total  energy  of  debris  ions  loss  down  field  lines 

-2  —1 

debris  ion  flux  (particles  cm  sec  ) 


COMMON/ I LEAK/ 
NLOSS : 

ELOSS: 

IFLUX  (2000): 


hour  of  HANE,  military  time 
latitude  of  HANE 

day  of  year  based  on  sequential  numbering 
altitude 

energy  yield  of  HANE  in  megatons 

weapon  mass  in  kilograms 

ionization  distance 

level  of  background  ionization 

total  number  of  debris  ions 

Radius  at  a  weapon  mass  of  ambient  plasma 


index  where  VIX  (1,2)  has  its  maximum 
index  where  VIX  (1,2)  has  its  minimum 
see  subroutine  rezone 
old  x-coordinates 


V.  CASE  RUN 


The  following  initial  data  and  boundary  condition  factors  were  used  to 
initialize  the  run  presented  below. 

The  output  for  this  run  is  presented  below.  No  physical  interpretation 
of  this  output  will  be  presented  here  since  the  results  are  only  presented  for 
code  documentation  purposes.  The  initialization  parameters  given,  while 
arbitrary,  are  designed  to  test  the  code  performance.  The  run  starts  at  time 
■  0.0  and  runs  up  to  0.5  millisecs.  Figure  1  illustrates  the  initial  profiles 
selected  (viz.,  densities,  temperatures,  pressures,  electric  fields,  and 
magnetic  field),  while  Figure  2  illustrates  the  formation  of  the  debris-air 
interface  t  a  0.5  millisecs.  Within  this  interface  air  is  swept  up  by  the 
debris,  a  laminar  electric  field  capable  of  debris  air  acceleration  is  formed, 
and  significant  ion-electron  heating  occurs. 


MINSTP  *  1. 

MAXSTP  *  101. 

I FLICK  *  500. 

ID  I AG  *  10, 

IPRNT  =  100, 

IPAGE  =  100, 

IDUMP  =  10, 

NSPEC  =  2, 

•END 
•BC I 

NIR  <  1  )  *  0.  0,  NIR  (  2  )  *  1.  0, 

NIL<  1  )  =1.0,  NIL <  2  )  *  1.  0, 

NIRO  <  1 )  =  5.  504E7,  NIR0(2>  =  0.0, 

N 1 LO  <  1 )  =  0.0,  NILO<  2 )  »  0.0, 

VIXR(l)  =  1  0,  VIXR<2>  *  1.0, 

VIXL<1)  =  10.  VIXL<2>  =  1.0, 

VI  XRO<  1  )  *  0  0.  VI  XRO  (  2  )  =  0.0, 

VIXL0<1)  =  00,  VI XLO ( 2 )  =  0.0, 

VI YR  <  1 )  =  l  .o,  VI YR  ( 2 )  =  1.0, 

VIYL(l)  *  1.0,  VIYL<2)  *  1.0, 

VIYRO(l)  =  00,  VIYR0<2)  *  0.0, 

VIYL0<1)  »  00.  VIYL0<2)  ■  0.0, 

VI2R<1)  =  10,  VI ZR  ( 2)  *  1.0, 

VI ZL  <  1 )  *  10.  VI ZL  ( 2 )  *  1.0, 

VIZRO(l)  =  0  0.  VI  ZRO  ( 2 )  *  0.0, 

VIZLO(l)  *  0.0,  VIZLCH2 )  ■  0.0, 

TIR  <  1  )  *0  0,  TIR  (2)  »  1.  0, 

TIL <  1)  «  i.  o,  TIL<2>  a  1.  0. 

TIRO ( l )  «  l.  161E7. TIRQ<2)  *  0.0, 

TILO(l)  =0  0.  TIL0(2>  *0.0, 

PIR(  1)  =0  0,  P IR  ( 2  )  *  1.  0, 

P IL <  1  )  *1.0,  P IL <2 )  »  1.  0. 

PIRO(l)  =  8  81E-2.  PIR0<2>  *  0.0, 

PILOID  *0.0,  PIL0<2>  a  o.  0, 

VIR  <  1 )  =0  0,  VIR  ( 2  )  a  0.  0, 

VIL  <  1 )  =0  0,  VIL  ( 2  )  *  0.  0. 

•END 
•  IOND 

MASS ( 1 )  =  2. 34E-23,  MASS<2>  =  4.  51E-23. 

ZZ  (  1  )  =10.  ZZ(2)  =1.0, 

TIMIN(l)  =  1  1605E4,  TIMIN<2)  =  00, 

NIMIN(l)  =  0  0.  NIMIN ( 2 >  *  0.0, 

TEMIN  *  1  1605E4, 

•END 

•BCE 


TER 

* 

0  0. 

TEL 

as 

1.  0, 

TERO 

a 

1.  161E6, 

TELO 

a 

0.  0. 

PER 

* 

1.  0, 

PEL 

* 

1.  0. 

PF.RO 

= 

0  0, 

PELO 

a 

0.  0. 

VER 

a 

0  0, 

VEL 

TZ 

0.  0, 

BYR 

a 

1. 0, 

BYL 

* 

1.  0. 

BYRO 

= 

0  0, 

BYLO 

* 

0.  0. 

BZR 

a 

0  0, 

BZL 

= 

1  0. 

3ZP0 

s 

0.  5, 

BZLO 

= 

0.  0. 

♦END 

•CONTRL 


LPAGE 

St 

F, 

LPRNT 

S 

T. 

LCPLE 

SB 

T, 

LCLASS 

LRESTR 

ac 

F, 

LDUMP 

SI 

T. 

LDIAG 

* 

T, 

LNEU 

LRAD 

s 

F, 

LANOM 

- 

T, 

LIAI 

s 

T. 

LBCI 

LMTI 

at 

T, 

LMI I 

X 

T. 

LUI I 

* 

T, 

LHDI 

LVISCO  -  T,  LDIFF  -  T,  LJOULE  *  F,  LARMOR  =  T, 

LSHUFF  *  T,  LI LEAK  =  F,  LTCENT  =  F,  LSTION  =  F, 

LDOHM  *  F,  LCHEM  *  T.  LJAC71  =  F,  LHANE  =  T, 

LSOLAR  *  F.  LLASER  *  F,  LECHO  =  T ,  LSTDP  =  F, 

LFLICK  »  F,  LREZON  *  T,  LCALCM  =  T, 

SEND 

♦INTDAT 

XMIN  *  1.  OE5, 

XMAX  ■  1.  0E6. 

NX  »  99, 

DT  =  1.  OE-IO. 

COURT  *  0.  2, 

ALPHA  -  3, 

TIMEO  *  0.  0, 

BYO  *  0  0, 

BYMIN  *  0.  O, 

BZO  *  0.  5, 

BZMIN  *  0  0, 

YSCALE  =  0  0, 

ZSCALE  =  0  0, 

VZERO  =  2.  0E8, 

NDZERO  -  1.  0E9, 

NAZERO  =  1.0E8, 

TDZERO  -  1.  1603E7. 

TAZERO  =  1.  1603E4, 

TEZERO  =  1.  1603E6. 

RSI  ZE  =  1.  OE 10, 

FLAG  =1.0, 

♦END 

♦HANE 

HRA  =  0.  0, 

LAT  =  40.  0. 

DAY  *  180, 

ALT  *  400  0. 

YIELD  =1.2, 

WMASS  =1.0, 

TEH  =  1.  0E6, 

BZRH  =  0.  5 

♦END 

♦SOLAR 

BETA  =  1.0E-3, 

VSOLAR  *  l .  OES, 

FMASS  =  1.0E12, 

FEWER  =  5.  0E27, 

TIAS  =  1.  0E6. 

TIDS  =  1. 0E4. 

TES  =  1.  0E6, 

NIAS  =  1.  OE 10, 

NIDS  =  1.  0E12, 

BZRS  =  1.  0E2, 

♦END 

♦LASER 

JOULES  =  30.  0, 

TMASS  *  4  51E-7. 

♦END 

♦CHEM 

NCHEM(l)  *  6, 

NCHEM  ( 2  >  =  9, 

ACHEM  =  1  0, 


♦END 
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HIKUUtS  CH-»»3 


(c) 

Figure  1.  Initial  conditions  (t  =  0).  Shown  are  the  (a)  densities,  (b) 
temperatures,  (c)  presures,  (d)  electric  fields,  and  (e)  magnetic 
field  profiles  used  to  initialize  the  run. 
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Figure  l.  Initial  conditions  (t  =0).  Shown  are  the  (a)  densities,  (b) 
(cont'd.)  temperatures,  (c)  pressures,  (d)  electric  fields,  and  (e) 
magnetic  field  profiles  used  to  initialize  the  run. 
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Figure  2.  Final  conditions  (t  =  5.3  x  10-^  sec).  Shown  are  the  (a) 

densities,  (b)  temperatures,  (c)  x  velocity,  (d)  y  velocity,  (e) 
pressures,  (f)  electric  fields,  and  (g)  magnetic  field  profiles 
after  approximately  0.5  millisecs. 
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DYNfcS  CM-*»2 


CX 


(g) 

Figure  2.  Final  conditions  (t  =  5.3  x  10""^  sec).  Shown  are  the  (a) 

(cont'd.)  densities,  (b)  temperatures,  (c)  x  velocity,  (d)  y  velocity,  (e) 

pressures,  (f)  electric  fields,  and  (g)  magnetic  field  profiles 
after  approximately  0.5  millisecs. 


25 


10000 1  tf-' 


Acknowledgments 


r 

I 


This  research  was  sponsored  by  the  Defense  Nuclear  Agency. 

References 

Boris,  J.P.:  1976,  "Flux  Corrected  Transport  Modules  for  Solving  Generalized 

Continuity  Equations,"  NRL  Memo  Rept.  3237.  (AD-A023-891) 

Braginskii,  S.I.:  1965,  "Reviews  of  Plasma  Physics  Vol.  I,”  p.  205,  Ed.  M. A. 

Leontovich,  Consultants  Bureau,  N.Y. 

Huba,  J.:  1983,  private  communication. 

Lampe,  M. ,  W.M.  Manheimer,  and  K.  Papadopoulos :  1975,  "Anomalous  Transport 

Coefficients  for  Hane  Applications  due  to  Plasma  Micro-instabilities,  NRL 
Memo  Rept.  3076.  (AD-A014-411) 


26 


DISTRIBUTION  LIST 


DEPARTMENT  OF  DEFENSE 

ASSISTANT  SECRETARY  OF  DEFENSE 
COMM,  CMD,  CONT  7  INTELL 
WASHINGTON,  D.C.  20301 

DIRECTOR 

COMMAND  CONTROL  TECHNICAL  CENTER 
PENTAGON  RM  BE  685 
WASHINGTON,  D.C.  20301 
01CY  ATTN  C-650 
01CY  ATTN  C-312  R.  MASON 

DIRECTOR 

DEFENSE  ADVANCED  RSCH  PROJ  AGENCY 
ARCHITECT  BUILDING 
1400  WILSON  BLVD. 

ARLINGTON,  VA.  22209 

01CY  ATTN  NUCLEAR  MONITORING  RESEARCH 
01CY  ATTN  STRATEGIC  TECH  OFFICE 

DEFENSE  COMMUNICATION  ENGINEER  CENTER 
1860  WIEHLE  AVENUE 
RESTON,  VA.  22090 
01 CY  ATTN  CODE  R410 
01CY  ATTN  CODE  R812 

DEFENSE  TECHNICAL  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA.  22314 
02CY 

DIRECTOR 

DEFENSE  NUCLEAR  AGENCY 
WASHINGTON,  D.C.  20305 
01CY  ATTN  STVL 
04CY  ATTN  TITL 
OICY  ATTN  DOST 
03CY  ATTN  RAAE 


DIRECTOR 

INTERSERVICE  NUCLEAR  WEAPONS  SCHOOL 
KIRTLAND  AFB,  NM  87115 

OICY  ATTN  DOCUMENT  CONTROL 

JOINT  CHIEFS  OF  STAFF 
WASHINGTON,  D.C.  20301 

OICY  ATTN  J-3  WWMCCS  EVALUATION  OFFICE 

DIRECTOR 

JOINT  STRAT  TGT  PLANNING  STAFF 
OFFUTT  AFB 
OMAHA,  NB  68113 
OICY  ATTN  JLTW-2 
OICY  ATTN  JPST  G.  GOETZ 

CHIEF 

LIVERMORE  DIVISION  FLD  COMMAND  DNA 
DEPARTMENT  OF  DEFENSE 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 
OICY  ATTN  FCPRL 

COMMANDANT 

NATO  SCHOOL  (SHAPE) 

APO  NEW  YORK  09172 

OICY  ATTN  U.S.  DOCUMENTS  OFFICER 

UNDER  SECY  OF  DEF  FOR  RSCH  &  ENGRG 
DEPARTMENT  OF  DEFENSE 
WASHINGTON,  D.C.  20301 

OICY  ATTN  STRATEGIC  &  SPACE  SYSTEMS  (OS) 

WWMCCS  SYSTEM  ENGINEERING  ORG 
WASHINGTON,  D.C.  20305 
OICY  ATTN  R.  CRAWFORD 

COMMANDER/DIRECTOR 
ATMOSPHERIC  SCIENCES  LABORATORY 
U.S.  ARMY  ELECTRONICS  COMMAND 
WHITE  SANDS  MISSILE  RANGE,  NM  88002 
OICY  ATTN  DELAS-EO  F.  NILES 


COMMANDER 
FIELD  COMMAND 
DEFENSE  NUCLEAR  AGENCY 
KIRTLAND,  AFB,  NM  87115 
OICY  ATTN  FCPR 

DEFENSE  NUCLEAR  AGENCY 
SAO/DNA 
BUILDING  20676 
KIRTLAND  AFB,  NM  87115 
OICY  D.C.  THORNBURG 


DIRECTOR 

BMD  ADVANCED  TECH  CTR 
HUNTSVILLE  OFFICE 
P.O.  BOX  1500 
HUNTSVILLE,  AL  35807 

OICY  ATTN  ATC-T  MELVIN  T.  CAPPS 
OICY  ATTN  ATC-0  W.  DAVIES 
01CY  ATTN  ATC-R  DON  RUSS 

PROCRAM  MANAGER 
BMD  PROGRAM  OFFICE 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 

01CY  ATTN  DACS-BMT  J.  SHEA 

CHIEF  C-E-  SERVICES  DIVISION 
U.S.  ARMY  COMMUNICATIONS  CMD 
PENTAGON  RM  1B269 
WASHINGTON,  D.C.  20310 

01CY  ATTN  C-  E-SERVICES  DIVISION 

COMMANDER 

FRADCOM  TECHNICAL  SUPPORT  ACTIVITY 
DEPARTMENT  OF  THE  ARMY 
FORT  MONMOUTH,  N.J.  .07703 

01CY  ATTN  DRSEL-NL-RD  H.  BENNET 
01CY  ATTN  DRSEL-PL-ENV  H.  BOMKE 
OICY  ATTN  J.E.  QUIGLEY 

COMMANDER 

U.S.  ARMY  COMM-ELEC  ENCRG  INSTAL  AGY 
FT.  HUACHUCA,  AZ  85613 

OICY  ATTN  CCC-EMEO  GEORGE  LANE 

COMMANDER 

U.S.  ARMY  FOREIGN  SCIENCE  &  TECH  CTR 
220  7TH  STREET,  NE 
CHARLOTTESVILLE,  VA  22901 
OICY  ATTN  DRXST-SD 

COMMANDER 

U.S.  ARMY  MATERIAL  DEV  &  READINESS  CMD 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 

OICY  ATTN  DRCLDC  J.A.  BENDER 

COMMANDER 

U.S.  ARMY  NUCLEAR  AND  CHEMICAL  AGENCY 
7500  BACKLICK  ROAD 
BLDG  2073 

SPRINGFIELD,  VA  22150 
OICY  ATTN  LIBRARY 


DIRECTOR 

U.S.  ARMY  BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN  PROVING  GROUND,  MD  21005 
OICY  ATTN  TECH  LIBRARY  EDWARD  BAICY 

COMMANDER 

U.S.  ARMY  SATCOM  AGENCY 
FT.  MONMOUTH,  NJ  07703 

OICY  ATTN  DOCUMENT  CONTROL 

COMMANDER 

U.S.  ARMY  MISSILE  INTELLIGENCE  AGENCY 
REDSTONE  ARSENAL,  AL  35809 
OICY  ATTN  JIM  GAMBLE 

DIRECTOR 

U.S.  ARMY  TRADOC  SYSTEMS  ANALYSIS  ACTIVITY 
WHITE  SANDS  MISSILE  RANGE,  NM  88002 
OICY  ATTN  ATAA-SA 
OICY  ATTN  TCC/F .  PAYAN  JR. 

OICY  ATTN  ATTA-TAC  LTC  J.  HESSE 

COMMANDER 

NAVAL  ELECTRONIC  SYSTEMS  COMMAND 
WASHINGTON,  D.C.  20360 

OICY  ATTN  NAVALEX  034  T.  HUGHES 
OICY  ATTN  PME  117 
OICY  ATTN  PME  117-T 
OICY  ATTN  CODE  5011 

COMMANDING  OFFICER 
NAVAL  INTELLIGENCE  SUPPORT  CTR 
4301  SUITLAND  ROAD,  BLDG.  5 
WASHINGTON,  D.C.  20390 

OICY  ATTN  MR.  DUBBIN  STIC  12 

OICY  ATTN  NISC-50 

OICY  ATTN  CODE  5404  J.  GALET 

COMMANDER 

NAVAL  OCCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152 
OICY  ATTN  J.  FERGUSON 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 


OICY 

ATTN 

CODE 

4700 

S. 

L.  Ossakow 

26  CYS  IF 

UNCLASS. 

1  CY  IF  CLASS) 

OICY 

ATTN 

CODE 

4701 

I 

Vitkov  It  sky 

OICY 

ATTN 

CODE 

4780 

J. 

Huba  (100 

CYS  IF  UNCLASS, 

1 

CY  IF  CLASS) 

OICY 

ATTN 

CODE 

7500 

OICY 

ATTN 

CODE 

7550 

OICY 

ATTN 

CODE 

7580 

OICY 

ATTN 

CODE 

7551 

OICY 

ATTN 

CODE 

7555 

OICY 

ATTN 

CODE 

4730 

E. 

MCLEAN 

OICY 

ATTN 

CODE 

4108 

OICY 

ATTN 

CODE 

4730 

B. 

RIPIN 

20CY 

ATTN 

CODE 

2628 

COMMANDER 

NAVAL  SEA  SYSTEMS  COMMAND 
WASHINGTON,  D.C.  20362 
01CY  ATTN  CAPT  R.  PITKIN 

COMMANDER 

NAVAL  SPACE  SURVEILLANCE  SYSTEM 
DAHLGREN,  VA  22448 

01CY  ATTN  CAPT  J.H.  BURTON 

OFFICER-IN-CHARGE 
NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK,  SILVER  SPRING,  MD  20910 
01CY  ATTN  CODE  F31 

DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20376 
OICY  ATTN  NSP-2141 
OICY  ATTN  NSSP-2722  FRED  WIMBERLY 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
DAHLGREN  LABORATORY 
DAHLGREN,  VA  22448 

OICY  ATTN  CODE  DF-14  R.  BUTLER 

OFFICER  OF  NAVAL  RESEARCH 
ARLINGTON,  VA  22217 
OICY  ATTN  CODE  465 
OICY  ATTN  CODE  461 
OICY  ATTN  CODE  402 
OICY  ATTN  CODE  420 
OICY  ATTN  CODE  421 


COMMANDER 

AEROSPACE  DEFENSE  COMMASD/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB ,  CO  80912 

OICY  ATTN  DC  MR.  LONG 

COMMANDER 

AEROSPACE  DEFENSE  COMMAND/XPD 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  80912 

OICY  ATTN  XPDQQ 

OICY  ATTN  XP 

AIR  FORCE  GEOPHYSICS  LABORATORY 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  OPR  HAROLD  GARDNER 
OICY  ATTN  LKB  KENNETH  S.W.  CHAMPION 
OICY  ATTN  OPR  ALVA  T.  STAIR 

OICY  ATTN  PHD  JURGEN  BUCHAU 

OICY  ATTN  PHD  JOHN  P.  MULLEN 

AF  WEAPONS  LABORATORY 
KIRTLAND  AFT,  NM  87117 
OICY  ATTN  SUL 

OICY  ATTN  CA  ARTHUR  H.  GUENTHER 
OICY  ATTN  NTYCE  1LT .  G.  KRAJEI 

AFTAC 

PATRICK  AFB,  FL  32925 
OICY  ATTN  TF/MAJ  WILEY 
OICY  ATTN  TN 

AIR  FORCE  AVIONICS  LABORATORY 
WRIGHT-PATTERSON  AFB,  OH  45433 
OICY  ATTN  AAD  WADE  HUNT 
OICY  ATTN  AAD  ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  &  ACQ 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  D.C.  20330 
OICY  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 
OICY  ATTN  J.  DEAS 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DI VISION/YSEA 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01732 
OICY  ATTN  YSEA 


29 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  DCKC  MAJ  J.C.  CLARK 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRICHT-PATTERSON  AFB,  OH  45433 
01CY  ATTN  NICD  LIBRARY 

01CY  ATTN  ETDP  B.  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
CRIFFISS  AFB,  NY  13441 

01CY  ATTN  DOC  LIBRARY/TSLD 
01CY  ATTN  OCSE  V.  COYNE 

SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

01CY  ATTN  SZJ 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 

OICY  ATTN  ADWATE  MAJ  BRUCE  BAUER 
01CY  ATTN  NRT 

OICY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.O.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

OICY  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CLAVTN 

SAMSO/MN 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

OICY  ATTN  MNNL 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB  MA  01731 

OICY  ATT:.  EEP  A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  C-042 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  FOR  A.  LABOWITZ 


DEPARTMENT  OF  ENERCY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

OICY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 

OICY  ATTN  DOC  CON  FOR  L-389  R.  OTT 

OICY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 

OICY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 


OICY 

ATTN 

DOC 

CON 

FOR 

J. 

WOLCOTT 

OICY 

ATTN 

DOC 

CON 

FOR 

R.I 

-.  TASCHEK 

OICY 

ATTN 

DOC 

CON 

FOR 

E. 

JONES 

OICY 

ATTN 

DOC 

CON 

FOR 

J. 

MALIK 

OICY 

ATTN 

DOC 

CON 

FOR 

R. 

JEFFRIES 

OICY 

ATTN 

DOC 

CON 

FOR 

J. 

ZINN 

OICY 

ATTN 

DOC 

CON- 

FOR 

P. 

KEATON 

OICY 

ATTN 

DOC 

CON 

FOR 

D. 

WESTER VELT 

OICY 

ATTN 

D.  ! 

SAPPENFIELD 

SANDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  W.  BROWN 

OICY  ATTN  DOC  CON  FOR  A.  THORNBROUCH 

OICY  ATTN  DOC  CON  FOR  T.  WRIGHT 

OICY  ATTN  DOC  CON  FOR  D.  DAHLGREN 

OICY  ATTN  DOC  CON  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT  DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  B.  MURPHEY 

OICY  ATTN  DOC  CON  FOR  T.  COOK 


OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERCY 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  DR.  YO  SONG 


OTHER  GOVERNMENT 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 

OICY  (ALL  CORRES:  ATTN  SEC  OFFICER  FOR) 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  4  INFO  ADMIN 
BOULDER,  CO  80303 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 

01CY  ATTN  W.  UTLAUT 
OICY  ATTN  D.  CROMBIE 
01CY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
OICY  ATTN  R.  GRUBB 
OICY  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  929S7 


LOS  ANGELES,  CA  90009 


OICY 

ATTN 

I.  CARFUNKEL 

OICY 

ATTN 

T.  SALMI 

OICY 

ATTN 

V.  JOSEPHSON 

OICY 

ATTN 

S.  BOWER 

OICY 

ATTN 

D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OICY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 
OICY  ATTN  L.  SLOAN 
OICY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 
OICY  ATTN  J.  WORKMAN 
OICY  ATTN  C.  PRETTIE 
OICY  ATTN  S.  BRECHT 


BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 
OICY  ATTN  G.  KEISTER 
OICY  ATTN  D.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
OICY  ATTN  D.B.  COX 
OICY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
OICY  ATTN  H.  LOCSTON 
OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  VJllshire  Blvd. 

Santa  Monica,  Calif  90401 
OICY  ATTN  C.B.  GABBARD 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
OICY  ATTN  F.  REIBERT 


31 


GENERAL  ELECTRIC  TECH  SERVICES  CO.,  ISC. 
HMES 

COURT  STREET 
SYRACUSE,  NY  13201 
01CY  ATTN  G.  MILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 

01 CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 

01CY  ATTN  TECHNICAL  LIBRARY 

01CY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01  CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
01CY  ATTN  J.M.  AEIN 
01CY  ATTN  ERNEST  BAUER 
01CY  ATTN  HANS  WOLFARD 
01CY  ATTN  JOEL  BENGSTON 

INTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

01CY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEMRA 
01CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01 CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
01CY  ATTN  DASIAC 
01CY  ATTN  WARREN  S.  KNAPP 
01CY  ATTN  WILLIAM  MCNAMARA 
01CY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
01CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
01CY  ATTN  DEPT  60-12 
01CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01CY  ATTN  MARTIN  WALT  DEPT  52-12 
01CY  ATTN  W.L.  IMHOF  DEPT  52-12 
01CY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
01CY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
01CY  ATTN  R.  HEFFNER 


M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

01CY  ATTN  DAVID  M.  TOWLE 
01CY  ATTN  L.  LOUGHLIN 
01CY  ATTN  D.  CLARK 
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MCDONNEL  DOUCLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CY  ATTN  N.  HARRIS 
OICY  ATTN  J.  HOULE 
01CY  ATTN  GEORGE  MROZ 
01CY  ATTN  W.  OLSON 
01CY  ATTN  R.W.  HALPRIN 
01CY  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01CY  ATTN  P.  FISCHER 
01CY  ATTN  W.F.  CREVIER 
OICY  ATTN  STEVEN  L.  GUTSCHE 
OICY  ATTN  R.  BOCUSCH 
OICY  ATTN  R.  HENDRICK 
OICY  ATTN  RALPH  KILB 
OICY  ATTN  DAVE  SOWLE 
OICY  ATTN  F.  FAJEN 
OICY  ATTN  M.  SCHEIBE 
OICY  ATTN  CONRAD  L.  LONGMIRE 
OICY  ATTN  B.  WHITE 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NEW  MEXICO  87106 
OICY  R.  STELLINGWERF 
OICY  M.  ALME 
OICY  L.  WRIGHT 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

OICY  ATTN  JOHN  MORGANSTERN 
OICY  ATTN  G.  HARDING 
OICY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
OICY  ATTN  W.  HALL 
OICY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

OICY  ATTN  E.C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  UA  98009 

OICY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  WILLIAM  B.  WRIGHT,  JR. 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  J.  KARZAS 
OICY  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 
OICY  ATTN  L.  DeRAND 
QICY  ATTN  W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICY  ATTN”  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

OICY  ATTN  VINCE  TRAPANI 


SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

01CY  ATTN  LEWIS  M.  LI.NSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  E.  FRIEMAN 
01CY  ATTN  E.A.  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  COODRIDCE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 


OICY 

ATTN 

DONALD  NETLSON 

OICY 

ATTN 

ALAN  BURNS 

OICY 

ATTN 

G.  SMITH 

OICY 

ATTN 

R.  TSUNODA 

OICY 

ATTN 

DAVID  A.  JOHNSON 

OICY 

ATTN 

WALTER  C.  CHESNUT 

OICY 

ATTN 

CHARLES  L.  RINO 

OICY 

ATTN 

WALTER  JAYE 

OICY 

ATTN 

J.  VICKREY 

OICY 

ATTN 

RAY  L.  LEADABRAND 

OICY 

ATTN 

G.  CARPENTER 

OICY 

ATTN 

G.  PRICE 

OICY 

ATTN 

R.  LIVINGSTON 

OICY 

ATTN 

V.  GONZALES 

OICY 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
01CY  ATTN  JOHN  ISE,  JR. 
01CY  ATTN  JOEL  GARBARINO 


TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
01CY  ATTN  R.  K.  PLEBUCH 
01CY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 
01CY  ATTN  D /  STOCKWELL 
SNTF/1575 

VIS1DYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
01CY  ATTN  W.  REIDY 
01CY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 
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